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Abstract 

In this talk we discuss mathematical structures associated to Feynman graphs. Feynman 
graphs are the backbone of calculations in perturbative quantum field theory. The mathemat- 
ical structures - apart from being of interest in their own right - allow to derive algorithms 
for the computation of these graphs. Topics covered are the relations of Feynman integrals 
to periods, shuffle algebras and multiple poly logarithms. 



1 Introduction 



High-energy particle physics has become a field where precision measurements have become 
possible. Of course, the increase in experimental precision has to be matched with more accurate 
calculations from the theoretical side. As theoretical calculations are done within perturbation 
theory, this implies the calculation of higher order corrections. This in turn relies to a large extent 
on our abilities to compute Feynman loop integrals. These loop calculations are complicated by 
the occurrence of ultraviolet and infrared singularities. Ultraviolet divergences are related to the 
high-energy behaviour of the integrand. Infrared divergences may occur if massless particles are 
present in the theory and are related to the low-energy or collinear behaviour of the integrand. 

Dimensional regularisation [1-3] is usually employed to regularise these singularities. Within 
dimensional regularisation one considers the loop integral in D space-time dimensions instead 
of the usual four space-time dimensions. The result is expanded as a Laurent series in the 
parameter £ = (4 — D)/2, describing the deviation of the D-dimensional space from the usual 
four-dimensional space. The singularities manifest themselves as poles in 1/e. Each loop can 
contribute a factor 1/e from the ultraviolet divergence and a factor 1/e 2 from the infrared di- 
vergences. Therefore an integral corresponding to a graph with / loops can have poles up to 
l/e 2/ . 

At the end of the day, all poles disappear: The poles related to ultraviolet divergences are 
absorbed into renormalisation constants. The poles related to infrared divergences cancel in the 
final result for infrared-safe observables, when summed over all degenerate states or are absorbed 
into universal parton distribution functions. The sum over all degenerate states involves a sum 
over contributions with different loop numbers and different numbers of external legs. 

However, intermediate results are in general a Laurent series in e and the task is to determine 
the coefficients of this Laurent series up to a certain order. At this point mathematics enters. We 
can use the algebraic structures associated to Feynman integrals to derive algorithms to calculate 
them. A few examples where the use of algebraic tools has been essential are the calculation of 
the three-loop Altarelli-Parisi splitting functions [4, 5] or the calculation of the two-loop ampli- 
tude for the process e + e~ — > 3 jets [6-15]. 

On the other hand is the mathematics encountered in these calculations of interest in its own 
right and has led in the last years to a fruitful interplay between mathematicians and physicists. 
Examples are the relation of Feynman integrals to mixed Hodge structures and motives, as well 
as the occurrence of certain transcendental constants in the result of a calculation [16-33]. 

This article is organised as follows: After a brief introduction into perturbation theory (sect. [2]), 
multi-loop integrals (sect. [3]> and periods (sect. U), we present in sect. |5] a theorem stating that 
under rather weak assumptions the coefficients of the Laurent series of any multi-loop integral 
are periods. The proof is sketched in sect.[6]and sect. [71 Shuffle algebras are discussed in sect.[8l 
Sect.|9]is devoted to multiple poly logarithms. In sect.[10]we discuss how multiple polylogarithms 
emerge in the calculation of Feynman integrals. Finally, sect.QTIcontains our conclusions. 
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2 Perturbation theory 



In high-energy physics experiments one is interested in scattering processes with two incoming 
particles and n outgoing particles. Such a process is described by a scattering amplitude, which 
can be calculated in perturbation theory. The amplitude has a perturbative expansion in the 
(small) coupling constant g: 

*n = g n (4 0) +8 2 4 l) +8 4 4 2) +8 6 4 3) + ..). (1) 

To the coefficient Sin contribute Feynman graphs with I loops and (n + 2) external legs. The 
recipe for the computation of Sin is as follows: Draw first all Feynman diagrams with the given 
number of external particles and / loops. Then translate each graph into a mathematical formula 

with the help of the Feynman rules. Sin is then given as the sum of all these terms. 

Feynman rules allow us to translate a Feynman graph into a mathematical formula. These 
rules are derived from the fundamental Lagrange density of the theory, but for our purposes 
it is sufficient to accept them as a starting point. The most important ingredients are internal 
propagators, vertices and external lines. For example, the rules for the propagators of a fermion 
or a massless gauge boson read 

n . . . // + m 

Fermion: = i- 



p 2 — m 2 + z'8 ' 



Gauge boson: WW = (2) 

k 2 + id 

Here p and k are the momenta of the fermion and the boson, respectively, m is the mass of the 
fermion. pf = p^ is a short-hand notation for the contraction of the momentum with the Dirac 
matrices. The metric tensor is denoted by g^ and the convention adopted here is to take the 
metric tensor as g^ v = diag(l, —1,-1,-1). The propagator would have a pole for p 2 = m 2 , or 
phrased differently for E = ±a/p 2 +m 2 . When integrating over E, the integration contour has 
to be deformed to avoid these two poles. Causality dictates into which directions the contour has 
to be deformed. The pole on the negative real axis is avoided by escaping into the lower complex 
half-plane, the pole at the positive real axis is avoided by a deformation into the upper complex 
half-plane. Feynman invented the trick to add a small imaginary part id to the denominator, 
which keeps track of the directions into which the contour has to be deformed. In the following 
we will usually suppress the j'8-term in order to keep the notation compact. 

As a typical example for an interaction vertex let us look at the vertex involving a fermion 
pair and a gauge boson: 



igf. (3) 
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Here, g is the coupling constant and y" denotes the Dirac matrices. At each vertex, we have 
momentum conservation: The sum of the incoming momenta equals the sum of the outgoing 
momenta. 

To each external line we have to associate a factor, which describes the polarisation of the 
corresponding particle: There is a polarisation vector £F{k) for each external gauge boson and a 
spinor u{p), u{p), v(p) or v(p) for each external fermion. 

Furthermore there are a few additional rules: First of all, there is an integration 



/ 



d 4 k 

(4) 



(2tc) 



for each loop. Secondly, each closed fermion loop gets an extra factor of ( — 1). Finally, each 
diagram gets multiplied by a symmetry factor 1/5, where S is the order of the permutation group 
of the internal lines and vertices leaving the diagram unchanged when the external lines are fixed. 

Having stated the Feynman rules, let us look at two examples: The first example is a scalar 
two-point one-loop integral with zero external momentum: 




k 



This integral diverges at k — > °° as well as at k — > 0. The former divergence is called ultraviolet 
divergence, the later is called infrared divergence. Any quantity, which is given by a divergent 
integral, is of course an ill-defined quantity. Therefore the first step is to make these integrals 
well-defined by introducing a regulator. There are several possibilities how this can be done, but 
the method of dimensional regularisation [1-3] has almost become a standard, as the calculations 
in this regularisation scheme turn out to be the simplest. Within dimensional regularisation one 
replaces the four-dimensional integral over the loop momentum by an D-dimensional integral, 
where D is now an additional parameter, which can be a non-integer or even a complex number. 
We consider the result of the integration as a function of D and we are interested in the behaviour 
of this function as D approaches 4. The original divergences will then show up as poles in the 
Laurent series in £ = (4 — D)/2. 

As a second example we consider a Feynman diagram contributing to the one-loop correc- 
tions for the process e + e~ — > qgq, shown in fig. ED At high energies we can ignore the masses of 
the electron and the light quarks. From the Feynman rules one obtains for this diagram (ignoring 
coupling and colour prefactors): 

-v(P4)fu( P5 )-L [ ^^u( Pl W(P2)^yMyMfv(P3). (6) 

Pm J y 2% ) k 2 Pn k i H 

Here, p n = Pi + P2, Pm =Pi+P2 + P3, h = h -p l2 , k 3 = k 2 -p3- Further $ (p 2 ) = ^ X {P2), 
where £ x (p2) is the polarisation vector of the outgoing gluon. All external momenta are assumed 
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Figure 1 : A one-loop Feynman diagram contributing to the process e + e — > qgq. 

to be massless: pf = for i = 1..5. We can reorganise this formula into a part, which depends 
on the loop integration and a part, which does not. The loop integral to be calculated reads: 

(21C)4*2*2*2' U) 

while the remainder, which is independent of the loop integration is given by 

-v( P4 )fu(p 5 )— — r u(pi)^(p2)p f i2yvypy^afv(p 3 ). (8) 

^123^12 

The loop integral in eq. © contains in the denominator three propagator factors and in the nu- 
merator two factors of the loop momentum. We call a loop integral, in which the loop momentum 
occurs also in the numerator a "tensor integral". A loop integral, in which the numerator is in- 
dependent of the loop momentum is called a "scalar integral". The scalar integral associated to 
eq. © reads 

(9) 

There is a general method [34, 35] which allows to reduce any tensor integral to a combination 
of scalar integrals at the expense of introducing higher powers of the propagators and shifted 
space-time dimensions. Therefore it is sufficient to focus on scalar integrals. Each integral can 
be specified by its topology, its value for the dimension D and a set of indices, denoting the 
powers of the propagators. 



3 Multi-loop integrals 

Let us now consider a generic scalar /-loop integral Iq in D = 2m — 2e dimensions with n prop- 
agators, corresponding to a graph G. For each internal line j the corresponding propagator in 
the integrand can be raised to a power V/. Therefore the integral will depend also on the num- 
bers Vi,...,v„. It is sufficient to consider only the case, where all exponents are natural numbers: 
Vj E N. We define the Feynman integral by 

ft r (v;) l , Dl n 

T tl ( 2^-lD/2 f'd D k r 1 
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with v = Vi + ... + V n . /u is an arbitrary scale, called the renormalisation scale. The momenta qj 
of the propagators are linear combinations of the external momenta and the loop momenta. The 
prefactors are chosen such that after Feynman parametrisation the Feynman integral has a simple 
form: 



v inn r I " „ i \ <7/ v -('+ 1 ) D / 2 

Xj>0 



The functions U and f depend on the Feynman parameters and can be derived from the topology 
of the corresponding Feynman graph G. Cutting / lines of a given connected /-loop graph such 
that it becomes a connected tree graph T defines a chord c(T,G) as being the set of lines not 
belonging to this tree. The Feynman parameters associated with each chord define a monomial 
of degree /. The set of all such trees (or 1-trees) is denoted by T\ . The 1-trees T E T\ define U as 
being the sum over all monomials corresponding to the chords C{T,G). Cutting one more line 
of a 1-tree leads to two disconnected trees (7i, T2), or a 2-tree. T2 is the set of all such pairs. The 
corresponding chords define monomials of degree l + Each 2-tree of a graph corresponds to a 
cut defined by cutting the lines which connected the two now disconnected trees in the original 
graph. The square of the sum of momenta through the cut lines of one of the two disconnected 
trees 7\ or T2 defines a Lorentz invariant 




*T = \ y Pi . (12) 



The function fo is the sum over all such monomials times minus the corresponding invariant. 
The function J is then given by 5o plus an additional piece involving the internal masses mj. In 
summary, the functions U and J are obtained from the graph as follows: 



u = e n *j 



j€C(T,G) 



H = 



n x j] (~ s ti) . 

(Ti,Ti)ei 2 u jec(T u G) 



7 = Fo+U Y^Xjm). (13) 
7=1 



4 Periods 

Periods are special numbers. Before we give the definition, let us start with some sets of numbers: 
The natural numbers N, the integer numbers Z, the rational numbers Q, the real numbers M 
and the complex numbers C are all well-known. More refined is already the set of algebraic 
numbers, denoted by Q. An algebraic number is a solution of a polynomial equation with rational 
coefficients: 

x n + a n -ix"- 1 + ••• + <*(> = 0, cij E Q. (14) 
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As all such solutions lie in C, the set of algebraic numbers Q is a sub-set of the complex numbers 
C. Numbers which are not algebraic are called transcendental. The sets N, Z, Q and Q are 
countable, whereas the sets R, C and the set of transcendental numbers are uncountable. 

Periods are a countable set of numbers, lying between Q and C. There are several equivalent 
definitions for periods. Kontsevich and Zagier gave the following definition [36]: A period is a 
complex number whose real and imaginary parts are values of absolutely convergent integrals of 
rational functions with rational coefficients, over domains in R" given by polynomial inequalities 
with rational coefficients. Domains defined by polynomial inequalities with rational coefficients 
are called semi-algebraic sets. 

We denote the set of periods by P. The algebraic numbers are contained in the set of periods: 
Q 6 P. In addition, P contains transcendental numbers, an example for such a number is 7i: 



The integral on the r.h.s. clearly shows that n is a period. On the other hand, it is conjectured that 
the basis of the natural logarithm e and Euler's constant Je are not periods. Although there are 
uncountably many numbers, which are not periods, only very recently an example for a number 
which is not a period has been found [37]. 

We need a few basic properties of periods: The set of periods P is a Q-algebra [36,38]. In 
particular the sum and the product of two periods are again periods. 

The defining integrals of periods have integrands, which are rational functions with rational 
coefficients. For our purposes this is too restrictive, as we will encounter logarithms as integrands 
as well. However any logarithm of a rational function with rational coefficients can be written as 



5 A theorem on Feynman integrals 

Let us consider a general scalar multi-loop integral as in eq. (fTTI) . Let m be an integer and set 
D = 2m — 2e. Then this integral has a Laurent series expansion in £ 




(15) 




(16) 



l G = £ CjE j . 



(17) 



Theorem 1: In the case where 



1. all kinematical invariants sj are zero or negative, 



2. all masses m, and p are zero or positive (p ^ 0), 



3. all ratios of invariants and masses are rational, 
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the coefficients Cj of the Laurent expansion are periods. 
In the special case were 

1 . the graph has no external lines or all invariants sj are zero, 

2. all internal masses mj are equal to fi, 

3. all propagators occur with power 1, i.e. v 7 = 1 for all j, 
the Feynman parameter integral reduces to 



xj>0 i=i 



)u~ D/2 (18) 



and only the polynomial U occurs in the integrand. In this case it has been shown by Belkale 
and Brosnan [39] that the coefficients of the Laurent expansion are periods. 

Using the method of sector decomposition we are able to prove the general case [40]. We 
will actually prove a stronger version of theorem 1 . Consider the following integral 

j = r d^dii-^x^fixr^mpM^. d9) 

x/> 1=1 v=i /;=i 

The integration is over the standard simplex. The a's, Z?'s, J's and /'s are integers. The P's are 
polynomials in the variables x\, x n with rational coefficients. The polynomials are required 
to be non-zero inside the integration region, but may vanish on the boundaries of the integration 
region. To fix the sign, let us agree that all polynomials are positive inside the integration region. 
The integral J has a Laurent expansion 

J = % c fi j - (2°) 

Theorem 2: The coefficients Cj of the Laurent expansion of the integral J are periods. 

Theorem 1 follows then from theorem 2 as the special case at = Vi — l,bi = 0,r = 2,P\ = U, 
p 2 = di +e/i = v - (/+ l)£>/2 and d 2 + ef 2 = ZD/2 - v. 

Proof of theorem 2: To prove the theorem we will give an algorithm which expresses each 
coefficient cj as a sum of absolutely convergent integrals over the unit hypercube with integrands, 
which are linear combinations of products of rational functions with logarithms of rational func- 
tions, all of them with rational coefficients. Let us denote this set of functions to which the 
integrands belong by M . The unit hypercube is clearly a semi- algebraic set. It is clear that 
absolutely convergent integrals over semi- algebraic sets with integrands from the set M are pe- 
riods. In addition, the sum of periods is again a period. Therefore it is sufficient to express each 
coefficient Cj as a finite sum of absolutely convergent integrals over the unit hypercube with inte- 
grands from M . To do so, we use iterated sector decomposition. This is a constructive method. 
Therefore we obtain as a side-effect a general purpose algorithm for the numerical evaluation of 
multi-loop integrals. 
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6 Sector decomposition 



In this section we review the algorithm for iterated sector decomposition [41-47]. The starting 
point is an integral of the form 



/ 



^5(i-£x ; )in^)n^w]" ; ' (2d 



xj>0 i=l v =1 / J =l 



where /Ji = ai + zbi and Xj = cj + edj. The integration is over the standard simplex. The a's, b's, 
c's and rf's are integers. The P's are polynomials in the variables x\, x n . The polynomials are 
required to be non-zero inside the integration region, but may vanish on the boundaries of the 
integration region. The algorithm consists of the following steps: 

Step 0: Convert all polynomials to homogeneous polynomials. 

Step 1: Decompose the integral into n primary sectors. 

Step 2: Decompose the sectors iteratively into sub-sectors until each of the polynomials is of 
the form 

P = x™K..x%*(c + P'(x)), (22) 

where c ^ and P'(x) is a polynomial in the variables Xj without a constant term. In this case the 
monomial prefactor jc'" 1 ...x%" can be factored out and the remainder contains a non-zero constant 
term. To convert P into the form (|22l) one chooses a subset S = {(Xi , a^} C {1 , ...n} according 
to a strategy discussed in the next section. One decomposes the ^-dimensional hypercube into k 
sub- sectors according to 

l k 1 k 

fd n x = H[d n xfl Q(x ai >x ai ). (23) 

In the Z-th sub-sector one makes for each element of S the substitution 

Xai = x ai x a . fori ^l. (24) 

This procedure is iterated, until all polynomials are of the form (|22)) . Fig.[2lillustrates this for the 
simple example S = {1,2}. Eq. (|23l) gives the decomposition into the two sectors x\ > x-i and 
X2 > x\. Eq. ((24)) transforms the triangles into squares. This transformation is one-to-one for all 
points except the origin. The origin is replaced by the line x\ = in the first sector and by the 
line X2 = in the second sector. Therefore the name "blow-up". 

Step 3: The singular behaviour of the integral depends now only on the factor 

f\x? +£hi . (25) 



(=1 
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X2 



•i'2 



Xl 



•i'l 



■r-i 



•i'l 



Figure 2: Illustration of sector decomposition and blow-up for a simple example. 



We Taylor expand in the integration variables and perform the trivial integrations 

l 



J dxx a+b£ 
o 



1 



a + l+be 



(26) 



leading to the explicit poles in 1/e. 



Step 4: All remaining integrals are now by construction finite. We can now expand all ex- 
pressions in a Laurent series in £ and truncate to the desired order. 

Step 5: It remains to compute the coefficients of the Laurent series. These coefficients con- 
tain finite integrals, which can be evaluated numerically by Monte Carlo integration. We imple- 
mentecfj] the algorithm into a computer program, which computes numerically the coefficients of 
the Laurent series of any multi-loop integral [45]. 



7 Hironaka's polyhedra game 

In step 2 of the algorithm we have an iteration. It is important to show that this iteration ter- 
minates and does not lead to an infinite loop. There are strategies for choosing the sub-sectors, 
which guarantee termination. These strategies [48-52] are closely related to Hironaka's polyhe- 
dra game. 

'The program can be obtained from http: //www.higgs .de/'stefanw/ software .html | 
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Figure 3: Illustration of Hironaka's polyhedra game. 

Hironaka's polyhedra game is played by two players, A and B. They are given a finite set M 
of points m = (mi, m n ) in N+, the first quadrant of W . We denote by A C R" the positive 
convex hull of the set M. It is given by the convex hull of the set 

\J(m + R n + ). (27) 

The two players compete in the following game: 

1. Player A chooses a non-empty subset S C {1, n}. 

2. Player B chooses one element i out of this subset S. 

Then, according to the choices of the players, the components of all (mi , m n ) G M are replaced 
by new points (m\, m' n ), given by: 

m'j = m h ifjV z \ 

m 'i = Y* m i~ c -> ^ 

where for the moment we set c = 1. This defines the set M'. One then sets M = M' and goes 
back to step 1 . Player A wins the game if, after a finite number of moves, the polyhedron A is of 
the form 

A = m + R n + , (29) 

i.e. generated by one point. If this never occurs, player B has won. The challenge of the polyhe- 
dra game is to show that player A always has a winning strategy, no matter how player B chooses 
his moves. A simple illustration of Hironaka's polyhedra game in two dimensions is given in 
fig- 121 Player A always chooses S = {1,2}. In ref. [45] we have shown that a winning strategy 
for Hironaka's polyhedra game translates directly into a strategy for choosing the sub-sectors 
which guarantees termination. 
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8 Shuffle algebras 



Before we continue the discussion of loop integrals, it is useful to discuss first shuffle algebras 
and generalisations thereof from an algebraic viewpoint. Consider a set of letters A. The set A is 
called the alphabet. A word is an ordered sequence of letters: 

w = hh...l k . (30) 

The word of length zero is denoted by e. Let K be a field and consider the vector space of words 
over K. A shuffle algebra Si on the vector space of words is defined by 

(hh---k) • (h+l—h) = J! ^o(l)^o(2)--^o(r)) (31) 

shuffles O 

where the sum runs over all permutations a, which preserve the relative order of 1, 2, ...,k and of 
k+ 1, ...,r. The name "shuffle algebra" is related to the analogy of shuffling cards: If a deck of 
cards is split into two parts and then shuffled, the relative order within the two individual parts is 
conserved. A shuffle algebra is also known under the name "mould symmetral" [53]. The empty 
word e is the unit in this algebra: 

e ■ w = w • e = w. (32) 
A recursive definition of the shuffle product is given by 

(hh-lk) ■ (h+l-lr) = H [(h-lk) ■ (h+l-lr)} +4+1 [(hh-h) ■ (lk+2-lr)] • (33) 

It is well known fact that the shuffle algebra is actually a (non-cocommutative) Hopf algebra [54]. 
In this context let us briefly review the definitions of a coalgebra, a bialgebra and a Hopf algebra, 
which are closely related: First note that the unit in an algebra can be viewed as a map from K 
to A and that the multiplication can be viewed as a map from the tensor product A ® A to A (e.g. 
one takes two elements from A, multiplies them and gets one element out). 

A coalgebra has instead of multiplication and unit the dual structures: a comultiplication A 
and a counit e. The counit is a map from A to K, whereas comultiplication is a map from A to A <g> 
A. Note that comultiplication and counit go in the reverse direction compared to multiplication 
and unit. We will always assume that the comultiplication is coassociative. The general form of 
the coproduct is 

A(a) = £a ; 0) ®af } , (34) 

where denotes an element of A appearing in the first slot of A <g> A and af^ correspond- 
ingly denotes an element of A appearing in the second slot. Sweedler's notation [55] consists in 
dropping the dummy index i and the summation symbol: 

A(a) = a (1) ®a (2) (35) 
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The sum is implicitly understood. This is similar to Einstein's summation convention, except 
that the dummy summation index i is also dropped. The superscripts (*) and ( 2 ) indicate that a 
sum is involved. 

A bialgebra is an algebra and a coalgebra at the same time, such that the two structures are 
compatible with each other. Using Sweedler's notation, the compatibility between the multipli- 
cation and comultiplication is expressed as 

A(a-b) = (a (1) -£ (1) )® (aW-b^Y (36) 

A Hopf algebra is a bialgebra with an additional map from A to A, called the antipode S , 
which fulfils 

a W -S (a (2) ) =S (a (1) ) -a (2) =0 forage. (37) 

With this background at hand we can now state the coproduct, the counit and the antipode 
for the shuffle algebra: The counit e is given by: 

e(e) = l, e(hl 2 ...l n )=0. (38) 

The coproduct A is given by: 

k 

A(hl 2 ...l k ) = £(/;+!•. 4)® (h...lj). (39) 
The antipode S is given by: 

s(hh...h) = (-1)* hh-i-hh. (40) 

The shuffle algebra is generated by the Lyndon words. If one introduces a lexicographic ordering 
on the letters of the alphabet A, a Lyndon word is defined by the property 

w<v (41) 

for any sub-words u and v such that w = uv. 

An important example for a shuffle algebra are iterated integrals. Let [a,b] be a segment of 
the real line and f\, fi, ... functions on this interval. Let us define the following iterated integrals: 

I(f h f 2 ,...,f k ;a,b) = J dtifi(h) J dt 2 f2(t 2 )... f dt k f k (t k ) (42) 

a a a 

For fixed a and b we have a shuffle algebra: 

I(flif2,-,f k ;a,b)-I(f k+h ...,f r ;a,b) = £ Kfa{l),fa{2), —,fa(r)'><*,b), (43) 

shuffles C7 
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Figure 4: Sketch of the proof for the shuffle product of two iterated integrals. The integral over 
the square is replaced by two integrals over the upper and lower triangle. 



where the sum runs over all permutations a, which preserve the relative order of 1,2, and 
of k+ 1, ...,r. The proof is sketched in fig. |4l The two outermost integrations are recursively 
replaced by integrations over the upper and lower triangle. 

We now consider generalisations of shuffle algebras. Assume that for the set of letters we 
have an additional operation 

(.,.) : A® A — >A, 

h®h^(luh), (44) 

which is commutative and associative. Then we can define a new product of words recursively 
through 

(hh-k) * (k+i-lr) = h [(h-k) * (k+i-lr)]+h+i [(hh-k) * (k+2-l r )] 

+(h,l k+ i)[(h-k)*(k+2-lr)}. (45) 

This product is a generalisation of the shuffle product and differs from the recursive definition of 
the shuffle product in eq. (|33l) through the extra term in the last line. This modified product is 
known under the names quasi-shuffle product [56], mixable shuffle product [57], stuffle product 
[58] or mould Symmetrel [53]. Quasi-shuffle algebras are Hopf algebras. Comultiplication and 
counit are defined as for the shuffle algebras. The counit e is given by: 

e(e) = l, e(hl 2 ...l n )=0. (46) 

The coproduct A is given by: 

k 

&(hh-k) = (47) 
7=0 

The antipode S is recursively defined through 

s(hh-k) = -hk...l k -Y,s(lj + i...l k )*(h...lj). (48) 

y=i 

An example for a quasi-shuffle algebra are nested sums. Let n a and ny be integers with n a < nt, 
and let f\, fi, ■■■ be functions defined on the integers. We consider the following nested sums: 

rib 'l-l 4-1-1 

S(fi,f2,...,f k ;n a ,n b ) = £ Mh) £ h{i 2 )... £ A(ijt) (49) 

i\=n a i2=>l a 'k="a 
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Figure 5: Sketch of the proof for the quasi-shuffle product of nested sums. The sum over the 
square is replaced by the sum over the three regions on the r.h.s. 

For fixed n a and n/, we have a quasi-shuffle algebra: 

S(fi,f2,—,fk',n a ,nb)*S(fk+i,...,f r ;n a ,nb) = 

£ /i( z 'i) s (h,---,fk;n a ,h - l)*S(f k+h ...,f r ;n a ,h - 1) 

h=n a 
tit, 

+ £ fk{jl) S(fhf2,-Jk;n a Ji-l)*S(fk + 2,-Jr;n a Ji-l) 

jl=n a 

+ £ /i(0/*(0 S(fz,..., fk) n a ,i- 1) *S(f k+2 ,...J r ;n a ,i- 1) (50) 

i=n a 

Note that the product of two letters corresponds to the point-wise product of the two functions: 

(f U fj)(n) = fi(n)fj(n). (51) 

The proof that nested sums obey the quasi-shuffle algebra is sketched in Fig. |5l The outermost 
sums of the nested sums on the l.h.s of (l50l) are split into the three regions indicated in Fig. |5] 

9 Multiple polylogarithms 

In the previous section we have seen that iterated integrals form a shuffle algebra, while nested 
sums form a quasi-shuffle algebra. In this context multiple polylogarithms form an interesting 
class of functions. They have a representation as iterated integrals as well as nested sums. There- 
fore multiple polylogarithms form a shuffle algebra as well as a quasi-shuffle algebra. The two 
algebra structures are independent. Let us start with the representation as nested sums. The 
multiple polylogarithms are defined by [59-62] 

i 1 >i 2 >-->ik>0 1 k 

The multiple polylogarithms are generalisations of the classical polylogarithms \A n (x), whose 
most prominent examples are 

Lh(x) = £ — = -ln(l -x), Li 2 (x) = £ X —, (53) 
11=1 11 h=i h 



15 



as well as Nielsen's generalised polylogarithms [63] 

Sn,p( X ) = L Wl,lr--,l( X '1^3)' ( 54 ) 

p-1 

and the harmonic polylogarithms [64, 65] 

^mi,...,mjfc(-*') = Li Wlr ^ Wfc (x, J, 1 ). (55) 

In addition, multiple polylogarithms have an integral representation. To discuss the integral 
representation it is convenient to introduce for Zk^O the following functions 



tk-i 

, r dn r dt 2 r dt k 
G{zh-,z k ;y) = / / ••• / 

J t\—Z\ J t2~Z2 J 



t\—Z\J t2~Z2 J tk — Zk 




(56) 



In this definition one variable is redundant due to the following scaling relation: 

G(zi,...,Zk;y) = G(xzi,...,xzk;xy) (57) 



If one further defines 



then one has 



g(z;y) = — , (58) 
y-z 



—G(zi,...,Zk;y) = g(zi;y)G(z2,---,z k ;y) (59) 
dy 



and 



G(zi,Z2,---,z k ;y) = f dt g(zi;t)G(z2,-,Zk;t). (60) 

o 

One can slightly enlarge the set and define G(0, ...,0;y) with k zeros for zi to Zk to be 

G(0,...,0;y) = ^(lny) k . (61) 

This permits us to allow trailing zeros in the sequence (zi,---,Zk) by defining the function G with 
trailing zeros via (1601) and (|6TT) . To relate the multiple polylogarithms to the functions G it is 
convenient to introduce the following short-hand notation: 

Gm u ...,m k (zi,—,Zk;y) = G(0^^,zi,...,z fc _i,a^0,z /t ;v) (62) 

mi— 1 mk~ 1 
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Here, all Zj for j = 1, ...,k are assumed to be non-zero. One then finds 



^-^m\,...,mk{x\i •••i^k) (1) ^m\,...,m\,_ ( ? >•••; »1 1 ■ (63) 

The inverse formula reads 

G mi mk (zi, ...,Zk',y) = (-l) fc Li raiv .. ;rajt (—,—,..., ) . (64) 

Eq. (|63l) together with (|62l) and (l56l) defines an integral representation for the multiple polyloga- 
rithms. 

Up to now we treated multiple polylogarithms from an algebraic point of view. Equally 
important are the analytical properties, which are needed for an efficient numerical evaluation. 
As an example I first discuss the numerical evaluation of the dilogarithm [66]: 

Li 2 (x) = - dt \ = £ - (65) 

J t n —\ n 



The power series expansion can be evaluated numerically, provided [jcj < 1 . Using the functional 
equations 

Li 2 (x) = - U 2(l)-^-\(H-x)f, 

U 2 (x) = -Li 2 (l-x) + — -ln(jc)ln(l-x). (66) 

6 

any argument of the dilogarithm can be mapped into the region |jc| < 1 and — 1 < Re(x) < 1/2. 
The numerical computation can be accelerated by using an expansion in [— ln(l —x)] and the 
Bernoulli numbers 5,-: 

Li 2W = IrTTTT(- ln (l-^))' +1 - (67) 

The generalisation to multiple polylogarithms proceeds along the same lines [67]: Using the 
integral representation eq. (1561) one transforms all arguments into a region, where one has a con- 
verging power series expansion. In this region eq. (|52l) may be used. However it is advantageous 
to speed up the convergence of the power series expansion. This is done as follows: The multiple 
polylogarithms satisfy the Holder convolution [58]. For z\ ^ 1 and z w ^ this identity reads 

G( Zl ,...,z w ;l)= (68) 
1 -zj,l — z>-i, 1 — zi; 1 -- ) g(zj+i,...,z w \- 



1 1 1 



P 



The Holder convolution can be used to accelerate the convergence for the series representation 
of the multiple polylogarithms. 
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10 From Feynman integrals to multiple polylogarithms 



In sect. |3] we saw that the Feynman parameter integrals depend on two graph polynomials Zl 
and J , which are homogeneous functions of the Feynman parameters. In this section we will 
discuss how multiple polylogarithms arise in the calculation of Feynman parameter integrals. 
We will discuss two approaches. In the first approach one uses a Mellin-Barnes transformation 
and sums up residues. This leads to the sum representation of multiple polylogarithms. In the 
second approach one first derives a differential equation for the Feynman parameter integral, 
which is then solved by an ansatz in terms of the iterated integral representation of multiple 
polylogarithms. 

Let us start with the first approach. Assume for the moment that the two graph polynomials 
11 and 7 are absent from the Feynman parameter integral. In this case we have 

JJ*^ )«<-£«> - r(vT + ... + v,r m 

With the help of the Mellin-Barnes transformation we now reduce the general case to eq. (|69l . 
The Mellin-Barnes transformation reads 



1 1 



Zoo loo 



(Ai+A 2 + ...+A n ) c = ^ - . / doi... I do n -i (70) 



- /CO 



Gn-i a -oi-...-a„_i-c 



xr(-c 1 )...r(-a n _i)r(c 1 + ... + a w _i+c) A?..A°rjAn 

Each contour is such that the poles of T(— a) are to the right and the poles of T(a + c) are to the 
left. This transformation can be used to convert the sum of monomials of the polynomials U and 
J into a product, such that all Feynman parameter integrals are of the form of eq. (|69l) . As this 
transformation converts sums into products it is the "inverse" of Feynman parametrisation. 

With the help of eq. (|69l) and eq. (1701 ) we may exchange the Feynman parameter integrals 
against multiple contour integrals. A single contour integral is of the form 

y+joo 

j = J_ r d r(a+ai)...r(o+fl w ) r(-c+fri)...r(-c+fr H ) 

y— ;'oo 

If max(Re(— ai), ...,Re(— a m )) < min(Re(&i), ...,Re(^„)) the contour can be chosen as a straight 
line parallel to the imaginary axis with 

max(Re(-ai),...,Re(-a m )) < Rey < min (Re(fci ),..., Re(fc„)) , (72) 

otherwise the contour is indented, such that the residues of T{p + a\), r(o + a m ) are to the 
right of the contour, whereas the residues of T(—a + bi), r(— o + b n ) are to the left of the 
contour. The integral eq. (1711) is most conveniently evaluated with the help of the residuum 
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theorem by closing the contour to the left or to the right. To sum up all residues which lie inside 
the contour it is useful to know the residues of the Gamma function: 

(-1)" (-1)" 
res (r(o + a),a= -a-n) = - — res (T(—o + a),a = a + n) = — - — (73) 

In general there are multiple contour integrals, and as a consequence one obtains multiple sums. 
Having collected all residues, one then expands the Gamma-functions: 

T(n + e)= (74) 
r(l+e)r(n) [l+£Z 1 (/i-l)+£ 2 Z 11 (n-l)+£ 3 Z 111 (n-l) + ...+£' , - 1 Z 11 ... 1 (n-l)] , 

where Z mir .. )Wjt (n) are Euler-Zagier sums defined by 

Zm u ...,m k {n) = J! Jm[---Jm^- ( 7 ^) 

n>h>i 2 >...>i k >0 1 k 

This motivates the following definition of a special form of nested sums, called Z-sums [68-71]: 

x h x h 

Z(n;m ll ...,m k ;x h ...,x k ) = £ ( 76 ) 

n>h>i 2 >...>i k >0 Zl lk 

k is called the depth of the Z-sum and w = m\ + . .. + is called the weight. If the sums go to 
infinity (n = °°) the Z-sums are multiple poly logarithms: 

Z(oo;mi,...,m k ;xi,...,x k ) = Li mu _^ mk (xi,...,x k ). (77) 

For x\ = ... = Xk = 1 the definition reduces to the Euler-Zagier sums [72-76]: 

Z(n;mi,...,m fe ;l,...,l) = Z mu _^ mk (n). (78) 

For n = o° and x\ = ... =xi c = 1 the sum is a multiple C,- value [58,77]: 

Z(oo;mi,...,TO*;l,...,l) = Cmi,...,mf (79) 

The usefulness of the Z-sums lies in the fact, that they interpolate between multiple polylog- 
arithms and Euler-Zagier sums. The Z-sums form a quasi-shuffle algebra. In this approach 
multiple poly logarithms appear through eq. (1771) . 

An alternative approach to the computation of Feynman parameter integrals is based on dif- 
ferential equations [65,78-82]. To evaluate these integrals within this approach one first finds 
for each master integral a differential equation, which this master integral has to satisfy. The 
derivative is taken with respect to an external scale, or a ratio of two scales. An example for a 
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one-loop four-point function is given by 




The two-point functions on the r.h.s are simpler and can be considered to be known. This equa- 
tion is solved iteratively by an ansatz for the solution as a Laurent expression in £. Each term in 
this Laurent series is a sum of terms, consisting of basis functions times some unknown (and to be 
determined) coefficients. This ansatz is inserted into the differential equation and the unknown 
coefficients are determined order by order from the differential equation. The basis functions are 
taken as a subset of multiple poly logarithms. In this approach the iterated integral representa- 
tion of multiple polylogarithms is the most convenient form. This is immediately clear from the 
simple formula for the derivative as in eq. (l59l) . 

11 Conclusions 

In this talk we reported on mathematical properties of Feynman integrals. We first showed that 
under rather weak assumptions all the coefficients of the Laurent expansion of a multi-loop inte- 
gral are periods. In the second part we focused on multiple polylogarithms and how they appear 
in the calculation of Feynman integrals. 

References 

[1] G. 't Hooft and M. J. G. Veltman, Nucl. Phys. B44, 189 (1972). 
[2] C. G. Bollini and J. J. Giambiagi, Nuovo Cim. B12, 20 (1972). 
[3] G. M. Cicuta and E. Montaldi, Nuovo Cim. Lett. 4, 329 (1972). 

[4] S. Moch, J. A. M. Vermaseren, and A. Vogt, Nucl. Phys. B688, 101 (2004), hep- 
ph/0403192. 



20 



[5] A. Vogt, S. Moch, and J. A. M. Vermaseren, Nucl. Phys. B691, 129 (2004), hep- 
ph/0404111. 

[6] L. W. Garland, T. Gehrmann, E. W. N. Glover, A. Koukoutsakis, and E. Remiddi, Nucl. 
Phys. B627, 107 (2002), hep-ph/01 12081. 

[7] L. W. Garland, T. Gehrmann, E. W. N. Glover, A. Koukoutsakis, and E. Remiddi, Nucl. 
Phys. B642, 227 (2002), hep-ph/0206067. 

[8] S. Moch, P. Uwer, and S. Weinzierl, Phys. Rev. D66, 114001 (2002), hep-ph/0207043. 

[9] A. Gehrmann-De Ridder, T. Gehrmann, E. W. N. Glover, and G. Heinrich, Phys. Rev. Lett. 

99, 132002 (2007), 0707.1285. 

[10] A. Gehrmann-De Ridder, T. Gehrmann, E. W. N. Glover, and G. Heinrich, JHEP 12, 094 
(2007), 0711.4711. 

[11] A. Gehrmann-De Ridder, T. Gehrmann, E. W. N. Glover, and G. Heinrich, Phys. Rev. Lett. 

100, 172001 (2008), 0802.0813. 

[12] A. Gehrmann-De Ridder, T. Gehrmann, E. W. N. Glover, and G. Heinrich, JHEP 05, 106 
(2009), 0903.4658. 

[13] S. Weinzierl, Phys. Rev. Lett. 101, 162001 (2008), 0807.3241. 

[14] S. Weinzierl, JHEP 06, 041 (2009), 0904.1077. 

[15] S. Weinzierl, Phys. Rev. D80, 094018 (2009), 0909.5056. 

[16] S. Bloch, H. Esnault, and D. Kreimer, Comm. Math. Phys. 267, 181 (2006), 
math. AG/05 10011. 

[17] S. Bloch and D. Kreimer, Commun. Num. Theor. Phys. 2, 637 (2008), 0804.4399. 

[18] S. Bloch, (2008), arXiv:0810.1313 [math.AG]. 

[19] F. Brown, Commun. Math. Phys. 287, 925 (2008), arXiv:0804.1660 [math.AG]. 

[20] F.Brown, (2009), arXiv:0910.01 14 [math.AG]. 

[21] F. Brown and K. Yeats, (2009), arXiv: 09 10.5429 [math-ph]. 

[22] O. Schnetz, (2008), arXiv:0801. 2856 [hep-th]. 

[23] O. Schnetz, (2009), arXiv: 0909. 0905 [math.CO]. 

[24] P. Aluffi and M. Marcolli, (2008), arXiv:0807.1690 [hep-th]. 

[25] P. Aluffi and M. Marcolli, (2008), arXiv:081 1.2514 [hep-th]. 

21 



[26] P. Aluffi and M. Marcolli, (2009), arXiv:0901.2107 [math.AG]. 

[27] P. Aluffi and M. Marcolli, (2009), arXiv:0907.3225 [math-ph]. 

[28] C. Bergbauer, R. Brunetti, and D. Kreimer, (2009), arXiv:0908.0633 [hep-th]. 

[29] S. Laporta, Phys. Lett. B549, 115 (2002), hep-ph/0210336. 

[30] S. Laporta and E. Remiddi, Nucl. Phys. B704, 349 (2005), hep-ph/0406160. 

[31] S. Laporta, Int. J. Mod. Phys. A23, 5007 (2008), 0803.1007. 

[32] D. H. Bailey, J. M. Borwein, D. Broadhurst, and M. L. Glasser, (2008), arXiv: 080 1.0891 
[hep-th]. 

[33] I. Bierenbaum and S. Weinzierl, Eur. Phys. J. C32, 67 (2003), hep-ph/03083 1 1. 
[34] O. V. Tarasov, Phys. Rev. D54, 6479 (1996), hep-th/9606018. 
[35] O. V. Tarasov, Nucl. Phys. B502, 455 (1997), hep-ph/9703319. 

[36] M. Kontsevich and D. Zagier, in: B. Engquis and W. Schmid, editors, Mathematics unlim- 
ited - 2001 and beyond , 771 (2001). 

[37] M. Yoshinaga, (2008), arXiv:0805.0349 [math.AG]. 

[38] B. Friedrich, (2005), arXiv:math.AG/0506113. 

[39] P. Belkale and P. Brosnan, Int. Math. Res. Not. , 2655 (2003). 

[40] C. Bogner and S. Weinzierl, J. Math. Phys. 50, 042302 (2009), 071 1.4863. 

[41] K. Hepp, Commun. Math. Phys. 2, 301 (1966). 

[42] M. Roth and A. Denner, Nucl. Phys. B479, 495 (1996), hep-ph/9605420. 

[43] T. Binoth and G. Heinrich, Nucl. Phys. B585, 741 (2000), hep-ph/0004013. 

[44] T. Binoth and G. Heinrich, Nucl. Phys. B680, 375 (2004), hep-ph/0305234. 

[45] C. Bogner and S. Weinzierl, Comput. Phys. Commun. 178, 596 (2008), 0709.4092. 

[46] A. V. Smirnov and M. N. Tentyukov, Comput. Phys. Commun. 180, 735 (2009), 0807.4129. 

[47] A. V. Smirnov and V. A. Smirnov, (2008), arXiv:08 12.4700. 

[48] H. Hironaka, Ann. Math. 79, 109 (1964). 

[49] M. Spivakovsky, Progr. Math. 36, 419 (1983). 

[50] S. Encinas and H. Hauser, Comment. Math. Helv. 77, 821 (2002). 

22 



[51] H. Hauser, Bull. Amer. Math. Soc. 40, 323 (2003). 
[52] D. Zeillinger, Enseign. Math. 52, 143 (2006). 

[53] J. Ecalle, (2002), (available at http://www.math.u-psud.fr/ biblio/ppo/2002/ppo2002- 
23.html). 

[54] C. Reutenauer, Free Lie Algebras (Clarendon Press, Oxford, 1993). 

[55] M. S weedier, Hopf Algebras (Benjamin, New York, 1969). 

[56] M. E. Hoffman, J. Algebraic Combin. 11, 49 (2000), math.QA/9907173. 

[57] L. Guo and W. Keigher, Adv. in Math. 150, 117 (2000), math.RA/0407155. 

[58] J. M. Borwein, D. M. Bradley, D. J. Broadhurst, and P. Lisonek, Trans. Amer. Math. Soc. 
353:3, 907 (2001), math.CA/9910045. 

[59] A. B. Goncharov, Math. Res. Lett. 5, 497 (1998), (available at http://www.math.uiuc.edu/K- 
theory/0297). 

[60] H. M. Minh, M. Petitot, and J. van der Hoeven, Discrete Math. 225:1-3, 217 (2000). 
[61] P. Cartier, Seminaire Bourbaki , 885 (2001). 

[62] G. Racinet, Publ. Math. Inst. Hautes Etudes Sci. 95, 185 (2002), math.QA/0202142. 
[63] N. Nielsen, Nova Acta Leopoldina (Halle) 90, 123 (1909). 

[64] E. Remiddi and J. A. M. Vermaseren, Int. J. Mod. Phys. A15, 725 (2000), hep-ph/9905237. 
[65] T. Gehrmann and E. Remiddi, Nucl. Phys. B601, 248 (2001), hep-ph/0008287. 
[66] G. 't Hooft and M. J. G. Veltman, Nucl. Phys. B153, 365 (1979). 

[67] J. Vollinga and S. Weinzierl, Comput. Phys. Commun. 167, 177 (2005), hep-ph/0410259. 

[68] S. Moch, P. Uwer, and S. Weinzierl, J. Math. Phys. 43, 3363 (2002), hep-ph/01 10083. 

[69] S. Weinzierl, Comput. Phys. Commun. 145, 357 (2002), math-ph/0201011. 

[70] S. Weinzierl, J. Math. Phys. 45, 2656 (2004), hep-ph/0402131. 

[71] S. Moch and P. Uwer, Comput. Phys. Commun. 174, 759 (2006), math-ph/0508008. 

[72] L. Euler, Novi Comm. Acad. Sci. Petropol. 20, 140 (1775). 

[73] D. Zagier, First European Congress of Mathematics, Vol. II, Birkhauser, Boston , 497 
(1994). 



23 



[74] J. A. M. Vermaseren, Int. J. Mod. Phys. A14, 2037 (1999), hep-ph/9806280. 

[75] J. Blumlein and S. Kurth, Phys. Rev. D60, 014018 (1999), hep-ph/98 10241. 

[76] J. Blumlein, Comput. Phys. Commun. 159, 19 (2004), hep-ph/031 1046. 

[77] J. Blumlein, D. J. Broadhurst, and J. A. M. Vermaseren, (2009), arXiv:0907.2557 [math- 
ph]. 

[78] A. V. Kotikov, Phys. Lett. B254, 158 (1991). 

[79] A. V. Kotikov, Phys. Lett. B267, 123 (1991). 

[80] E. Remiddi, Nuovo Cim. A110, 1435 (1997), hep-th/971 1 188. 

[81] T. Gehrmann andE. Remiddi, Nucl. Phys. B580, 485 (2000), hep-ph/99 12329. 

[82] T. Gehrmann and E. Remiddi, Nucl. Phys. B601, 287 (2001), hep-ph/0101 124. 



24 



